Holonomy for Quantum Channels 
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A quantum holonomy reflects the curvature of some underlying structure of quantum mechanical 
systems, such as that associated with quantum states. Here, we extend the notion of holonomy to 
families of quantum channels, i.e., trace-preserving completely positive maps. By the use of the 
Jamiolkowski isomorphism, we show that the proposed channel holonomy is related to the Uhlmann 
holonomy. The general theory is illustrated for specific examples. We put forward a physical 
realization of the channel holonomy in terms of interferometry. This enables us to identify a gauge- 
invariant physical object that directly relates to the channel holonomy. Parallel transport condition 
and concomitant gauge structure are delineated in the case of smoothly parametrized families of 
channels. Finally, we point out that interferometer tests that have been carried out in the past 
to confirm the An rotation symmetry of the neutron spin, can be viewed as early experimental 
realizations of the channel holonomy. 
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I. INTRODUCTION 

Quantum geometric phases and quantum holonomies 
have, since their initial discovery, proven to be a versa- 
tile structure that appears in many different contexts in 
quantum mechanics. Berry showed in his seminal paper 
[l| that a state vector initially in an eigenspace of a non- 
degenerate Hamiltonian acquires a geometric phase fac- 
tor in addition to the familiar dynamical phase factor af- 
ter being adiabatically transported around a closed curve 
in the parameter space of the Hamiltonian. Wilczek and 
Zee soon thereafter extended Berry's work by show- 
ing that the geometric phase factor generalizes to a uni- 
tary state change, called a non-Abelian holonomy (or just 
holonomy for short), in the case of degenerate Hamilto- 
nians. Another extension of Berry's work was provided 
by Aharonov and Anandan [3] who removed the require- 
ment of adiabaticity by showing that a geometric phase 
factor is defined for any cyclic evolution of a pure quan- 
tum state. This result was further generalized by Samuel 
and Bhandari [3| to include noncyclic evolution as well. 
A holonomy for curves of density operators was first in- 
troduced by Uhlmann 5]. 

The aforementioned geometric phases and holonomies 
may be classified in the following way. Firstly, we have 
holonomies for subspaces, such as eigenspaces of Hamil- 
tonians l|, subspaces selected by projective measure- 
ments m , cyclic subspaces 0] , and decoherence free sub- 
spaces H, . Holonomic quantum computation [13, [HI is 
related to this class of holonomies. Secondly, we have ge- 
ometric phases and holonomies for quantum states, both 
pure 0,13 and mixed |5[. The geometric structures of 
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these two classes are given by the fiber bundles associ- 
ated with the mappings, "basis of subspace" ^-^ "sub- 
space" and "purification" ^ "state", respectively. 

In this paper, we focus on the geometry related to a 
third major concept in quantum theory, namely quantum 
maps. More precisely, we are interested in the holonomy 
for sequences of quantum channels, i.e., trace-preserving 
completely positive maps. This concept, which we shall 
call "channel holonomy" , is associated with the geometry 
given by a fiber bundle structure related to the mapping 
"Kraus representation" i-^- "completely positive map" . 
The aim with this analysis is to delineate this structure 
and to examine its physical relevance. 

The outline of this paper is as follows. In the next 
section, we introduce the channel holonomy and exam- 
ine its behavior under gauge transformations. The rela- 
tion between the Uhlmann holonomy for a sequence 
of density operators, constructed from the Jamiolkowski 
isomorphism [T^j , and the channel holonomy is analyzed 
in Sec. IIIII The case of smoothly parametrized families of 
quantum channels is discussed in Sec. IIVI We derive the 
parallel transport condition and introduce a gauge poten- 
tial associated with such families. In Sec. |Vl the channel 
holonomy is calculated for specific types of channel se- 
quences. A physical realization of the channel holonomy 
based on ancillary constructions in two-beam interfer- 
ometry is demonstrated in Sec. IVII In particular, we 
demonstrate that the channel holonomy is related to the 
"gluing matrix" [l^, [iJl that arises when two channels 
are combined in an interferometer. The interferometer 
setup also provides means to identify a physically mean- 
ingful gauge-invariant object associated with the channel 
holonomy. In Sec. IVIH we examine the case of smoothly 
parametrized families of quantum channels using the an- 
cillary construction. It should be noted that the analysis 
in Sees. IVII and lVIII parallel to a large extent that of Sees. 
Ulland HVl the main difference being that while the latter 
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utilizes directly the Kraus operators, the former utilizes 
ancillary systems. In Sec. [ YHIi we examine the inter- 
ferometer tests in Refs. [1^, of the 47r symmetry of 
the neutron spin in terms of the channel holonomy. The 
paper ends with the conclusions. 

II. CHANNEL HOLONOMY 

Consider a trace-preserving completely positive map 
(channel for short) T acting on a Z?-dimensional state 
space of a quantum system. The action of the channel 
on a state p can be expressed as 

K 

Hp) = Y^P^pPl (1) 
fc=i 

The operators Fk constitute a Kraus representation of 
T We assume that K is the number of linearly 

independent Kraus operators needed to represent i.e., 
K is the Kraus number K{T) of the channel [l^. From 
trace preservation it follows that ^^.F^Fk = 1, where 1 
is the identity operator. The Kraus representation of a 
channel is not unique. If {Fk}k is a valid representation 
of !F then so is {Fk}k, where 

K 

Ffc = ^Fii7ife, (2) 

1=1 

U being a unitary matrix p^ . 

Let {Ek}k be a linearly independent Kraus represen- 
tation of a channel £. Given another channel J- with 
K{J^) — K{£) — K, we wish to find a linearly indepen- 
dent Kraus representation {Fk}k of J-^ that in some sense 
is parallel with {Ek}k- A convenient choice would be to 
find the Kraus representation {Fk}k that minimizes [l9l | 

J2\\Ek- Fkf = 2D- 2ReTrT, (3) 

k 

where || • j] denotes the Hilbert-Schmidt norm and T is a 
matrix with elements Tki = Ti{fIEi). Under a change 
of Kraus representation {Fk}k {Fk}k, as given by Eq. 
([21), the matrix T transforms as T -> T U^T. Hence, 

J2\\Ek-Fkf ^2D-2ReTT{U^T). (4) 

fe 

Assuming that T is of rank K (see the appendix for 
an elaboration on the rank of T) the minimum is ob- 
tained when T > which corresponds to the choice 

U ^ <^iT) = VtT^ T. Thus, we say that two Kraus 
representations {Ek}k and {Fk}k are parallel if their cor- 
responding T matrix is positive definite. 

We are now in a position where we can define a holon- 
omy corresponding to a sequence £i, . . . ,£n oi quantum 
channels with K{£i) = . . . = K{£m) = K . This is done 



by choosing a Kraus representation {E]^}^^^ for each £"„ 
in the sequence and encoding the Kraus freedom in a 
family of unitary matrices [/„, n = l,...,iV. In this 
way, we may express the parallelity conditions as 

C/|;+iT„+i,„i7„ > 0, n = l,...,iV-l, (5) 

where the K x K matrices Tn+i,n with elements 

[Tn+i.n]ki = TriE^^+'^E'n (6) 

are all assumed to be of rank K. The conditions in Eq. 
([5]) are satisfied by 

Un+i = $(T„+i,„)(7„. (7) 

We obtain after iteration 

Un ^HTn.n-i)---HT2^i)Ui. (8) 

Define 

U,i,{£i,. . . ,£n) = HTi^n)UnU\ 

= $(Ti,jv)$(TAr,Ar_i)...$(r2,l) (9) 

to be the channel holonomy for the sequence £i^ . . . ,£n- 
If we consider the set of linearly independent Kraus 
representations as a fiber bundle with the set of chan- 
nels with a fixed Kraus number as base manifold, the 
change of Kraus representations, as in Eq. can be 
interpreted as a gauge transformation. As seen from the 
definition in Eq. © the factor UnU\ is multiplied from 
the left by ^[Ti^n]. This construction guarantees that 
the matrix Uch[£i, ■ ■ ■ , £n] transforms gauge covariantly. 
To see this we make the gauge transformation 

K->E^r[v„];fc (10) 
I 

of the Kraus representation of each channel in the se- 
quence. We obtain 

$(T„+i,„) yt^i$(r„+i,„)y„, (ii) 

which implies 

i7ch(f 1, ■■■,£n)^ V\u,h{£u • • ■ , £n)Vi. (12) 

Hence, the channel holonomy transforms gauge covari- 
antly as required. 

Concerning the gauge covariance, let us point out that 
the channel holonomy is described as a unitary matrix. 
The gauge covariance of this matrix is necessary since 
it, in some sense, is a matrix representation of a gauge- 
invariant object, and the gauge covariance reflects the 
freedom of the choice of "basis" in this matrix represen- 
tation. It is to be noted [2^ that this is the case also for 
other type of holonomies, such as those of Refs. 0, [2]| . 
In Sec. IVII we elucidate what gauge-invariant object the 
channel holonomy represents, and in what sense the chan- 
nel holonomy is a matrix representation of this object. 
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An alternative way to obtain the channel holonomy 
defined above is to make a gauge transformation yielding 
Tn+i,n T'n+i,n > 0, for all n = 1,...,7V - 1, i.e., 
^{Tn+i,n) —1,1 being the K x K identity matrix. This 
choice amounts to parallel transport along the sequence 
£i, . . . ,fjv. From Eq. ^ it follows that 

[/ch(fl,...,fAr) = $(Ti,^), (13) 

which is the expression for the channel holonomy in the 
parallel transport gauge. 



III. RELATION TO THE UHLMANN 
HOLONOMY 

For any given channel £ acting on elements of £(7ig), 
i.e., the set of linear operators on the £)-dimensional 
Hilbert space Tig, one can find a corresponding den- 
sity operator p g ^{Ti-q ® Tiq). This can be done 
via the Jamiolkowski isomorphism [12| . i.e., £ ^ p — 
£ ®T{\'ip){il)\), where I is the identity channel and IV') = 
X^fcLi 1^) ® \k) with {|fc)}fc an orthonormal basis of 
Tiq. We show that the holonomy associated with a se- 
quence of channels f i , . . . , £iq is related to the Uhlmann 
holonomy 5] for the extended sequence of density oper- 
ators pi, . . . ,pn,Pn+i = Pi 22], where 

Pn^£n®im{^\)- (14) 

To avoid technical complications, we assume that all 
channels have maximal Kraus number K(£n) — . This 
guarantees that the corresponding density operators p„ 
are faithful 0], i.e., they are full rank. 

There are several ways to calculate the Uhlmann holon- 
omy associated with a sequence pi, . . . , pn+i of faithful 
density operators. For each density operator p„ there 
corresponds Uhlmann amplitudes Wn — \fp7iVn, where 

Vn = <&(W„) = \jWnWX Wn is a unitary operator. A 

sequence of such amplitudes is parallel if W\j^-^Wn > 0, 
which allows us to define the Uhlmann holonomy J/uhi — 
Vn+iVi . If we instead consider another sequence of am- 
plitudes Wn that is not parallel transported (but corre- 
sponds to the same sequence of density operators), we 
can make it into a parallel transported sequence WnUn 
by a choice of unitary operators C/„ such that 

uU^Wi^^WnUn > 0, (15) 

which implies that Un+i = ^{Wn^iWn)Un and we find 
Uuu = ^WN+i)UN+iUl<i>iWl). If we furthermore 
assume that the sequence of amplitudes is cyclic, i.e., 
Wn+1 = Wi (which implies that the underlying sequence 
of density operators is cyclic, i.e., pn+i = Pi), we obtain 

C/uhi = ^Wi)UN+iUl<i>{wl). (16) 



We also note that 

Un+1 = ^W^WmMWI^Wn-i) ■ ■ ■ HW^Wi)Ui. (17) 

Let us now return to the cyclic sequence of density 
operators pi, . . . , pjv, Pn+i = Pi as defined by Eq. (|14p . 
We fix an orthonormal basis {\fk)}k of Hq (8) Hq, and 
define the amplitudes 

Wn = J2(Ek®m){fk\, (18) 

k 

where {EJ^}k is a linearly independent Kraus representa- 
tion of £n ■ We furthermore find that 

wUiWn = ;^El/fc)(/'|Tr(i?ri^£;r) 

kl 

It follows that the two amplitudes Wn and Wn+i defined 
in Eq. (fT8|) are parallel if and only if Tn+i,n > 0. Hence, 
the parallelity condition for channels is, via this construc- 
tion, closely related to the Uhlmann parallel transport 
of amplitudes. If we combine Eq. (fT7|) with the fact 

that HWl^.Wn) = Eki\fk){fimTn+i,n)]ki wc Can 
conclude that 

f7w+i=^|/fc)(//|[C/ch]MC/i, (20) 

kl 

which implies 

[t/ch]fe; = {fkmwl)UvhMWi)\fi). (21) 

The right-hand side of the above equation transforms as 
in Eq. p2)) under the gauge transformation in Eq. (jlOp . 
as it should. We note that it is <I>(Vl^i) and ^{W^) that 
are "responsible" for the gauge covariance, as J7uhi re- 
mains invariant. Note also that the construction in Eq. 
(fT8)) leading to Eq. pT|) contains an arbitrary choice of 
basis {\fk)}k, as well as an arbitrary maximally entan- 
gled state The channel holonomy Uch should not 
depend on any of these arbitrary choices. (This invari- 
ance of Uch with respect to the choice of {\fk)}k and 
should not be confused with the covariance of Uch under 
the gauge transformations in Eq. PH]).) The Uhlmann 
holonomy J7uhi is invariant under change of {\fk)}k, and 
although $(Wi) depends on {|/fc)}fc, the combination 
<f>(Wi)|/fe) does not, which leaves Uch invariant under 
the choice of Given a maximally entangled state 

IV'), all other maximally entangled states can be obtained 
as [1 (g) U]\ip), where U is unitary. Both t/yhi and <I>(Wi) 
depend on the choice of \tjj). However, the combination 
^{Wl)Uuh\^{Wi) is independent of the choice of max- 
imally entangled state \tp). We can thus conclude that 
Uch behaves as required. 
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IV. SMOOTHLY PARAMETRIZED FAMILIES 
OF CHANNELS 

A. The parallel transport condition 

So far we have considered the holonomies associated 
with sequences of channels. Here we consider the tran- 
sition to families of channels £s smoothly parameter- 
ized by a real variable s e [0,1]. As before we as- 
sume that the Kraus number is constant within each 
family, i.e., K{Ss) = K for all s. Consider a smoothly 
parametrized family of linearly independent Kraus repre- 
sentations {Ek{s)}^^i of Eg. For sequences of Kraus rep- 
resentations, the parallel transport condition is Tn+i,n > 
0. A generalization of this condition to smooth curves 
is obtained by requiring that the matrix with elements 
Tt[eI{s + Ss)Ei{s)] is positive definite to first order in 
the limit of small Ss. We find 

Tt[eI{s + Ss)Eiis)] = + SsRki, (22) 

where 

Qki = Tt[eI{s)Ei{s)], Rki = Tt[eI{s)Ei{s)]. (23) 

Since {Ek{s)}k is a linearly independent set it follows 
that Q is positive definite. Therefore, a necessary con- 
dition for Q + 5sR to be positive definite is that R is 
Hermitian. We shall now see that this is also a suffi- 
cient condition for small 6s. Since Q is positive definite 
it follows that 

Q + SsR= ./Q{I + Ss^~^Ry/Q'^)y/Q. (24) 

The assumption R^ = R implies that ^/Q ^ R^/Q ^ is 
a finite Hermitian matrix and consequently its eigenval- 
ues are real and Q -I- 5sR is positive definite for suffi- 
ciently small 5s. We can conclude that a smooth fam- 
ily of Kraus representations {Ek{s)}^^^ of channels with 
constant Kraus number K is parallel transported if and 
only if 

i:T[El{s)Ei{s)]=i:r[El{s)Ei{s)], k,l^l,...,K. 

(25) 

This can be compared with the Uhlmann parallel trans- 
port condition for smoothly parameterized families of 
amplitudes W{s), which is W''{s)W{s) = W^{s)W{s) 
i- 

B. The gauge potential 

We now demonstrate how to introduce a gauge po- 
tential along the path C : [0, 1] 9 s i-^ £s of smoothly 
parametrized channels. Suppose we choose a family of 
Kraus representations {i?fc(s)}^]^ over C that is not par- 
allel transported. We can make a gauge transformation 
of the form 

K 

£;fe(s) ^^fc(s) = ^S/(s)!7,fc(s) (26) 
1=1 



such that {Ek{s)}^^i is parallel transported. By insert- 
ing Eq. (gHI) into Eq. we obtain 

R R'^^UU^Q + QUU\ (27) 

where Q and R arc as in Eq. ((23l) . Let A{s) be an anti- 
Hermitian matrix generating U{s) via U = AU, and 
substitute into Eq. ^ yielding (cf. Rcf. [23]) 

R-R'' = AQ + QA. (28) 

If we assume Q{s) > for all s, then, according to The- 
orem VH.2.3 in Ref. [li], we find that Eq. ^ has a 
unique solution A(s) for each s, namely 

/•OO 

A{s)= e-'-Q^''>[R{s)-R^s)]e-^^^'Ur. (29) 
Jo 

One can see from the right-hand side of this equation that 
A{s) is indeed an anti-Hermitian matrix for all s S [0, 1]. 

To see that A transforms as a proper gauge poten- 
tial [25j we consider an arbitrary gauge transformation 
Ek{s) ^-> i?;(s)y;fc(s), where V{s) is a smooth fam- 
ily of unitary operators. We find that Q V^QV and 

R V^QV + V^RV. From these transformation prop- 
erties we obtain 

A^V^AV + vW, (30) 
i.e., A indeed transforms as a gauge potential. 

C. An example 

The gauge potential A{s) along C can be found by 
solving Eq. (pS)) . Although this is difficult in the general 
case, it can be done for Kraus number K = 2. Let a = 
{(Tx, (Ty,(Tz) be the standard 2x2 Pauli matrices. Since 
Q > and TrQ = Tr/, we may write 

Q = I + x a, (31) 

where ]x] ^ 1 (the eigenvalues of Q > are 1 ± ]x]). 
Furthermore, from the definition of the matrices Q and 
R, it follows that TrR + TrR^ = TrQ = 0, which implies 
RcTrJ? = 0. Thus, we may put 

R = izol + (y + iz) ■ a, (32) 

where 2y = x. Finally, A is anti-Hcrmitian, which sug- 
gests that we can write 

A = iuol + iu ■ (T. (33) 

We assume that x, zo,y,z,uo,u are all smooth real- 
valued functions of s. By inserting Eqs. (|3ip. (15^ . and 
((33| into Eq. we obtain 



/ Z — ZnX -I- X X (x X z) \ , , 

( 1 - Ix]^ -)-'^ (^4) 
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along C, which is well-defined since |x| ^ 1. Note, in 
particular, that A vanishes when and z vanish, which 
from Eq. ([32| can be seen to correspond to the parallel 
transport condition i?^ = R. Thus, the form of the K = 
2 gauge potential is consistent with the notion of parallel 
transport developed in Sec. IIV Al Note also that a result 
analogous to Eq. for the Uhlmann holonomy in the 
case of smooth families of faithful qubit density operators 
has been obtained in Ref. 1261. 



D. Writing the holonomy in terms of the gauge 
potential 

It remains to demonstrate that the holonomy of the 
smooth path C : [0, 1] 9 s i— s- fj. can be expressed in 
terms of A. Consider the polar decomposition Tg^gg s — 
\T s+Ss,s\^{T s+Ss,s) ■ To the first order in Ss we may write 
|T5+5s,s| = Q + SsH and <^{Ts+Ss,s) = I + SsJ, where 
H'^ = H and = - J. Thus, T,+s,^s = Q + Ss{QJ + 
H) + 0{Ss^) and we obtain 

T,+ss,.-TI^s,^^ = ds{QJ + JQ) + 0i5s^). (35) 



Furthermore, from Eq. ([28]) and the expression Ts+Ss,s — 
Q + SsR + OiSs^), we &nd 

Ts+ss,s - Tl^ss.^ ^ Ss{QA + AQ) + OiSs^). (36) 

Since Q > it follows that J — A and the channel 
holonomy can be written as 

i7ch(C) =$(To,i)Pe^o (37) 
where P denotes path ordering and ^ is a solution of Eq. 



V. EXAMPLES 

In this section, we consider two examples where the 
channel holonomy can be explicitly calculated. The first 
example concerns unitary channels, i.e., operations on 
closed quantum systems. It turns out that the chan- 
nel holonomy for such channels is intimately connected 
to certjiin cases of the mixed state phase proposed in 
Ref. [27j. Our second example concerns channels for 
which there exist Kraus representations that are built up 
by subspace holonomies associated with smooth paths 
that, e.g., can be approximately generated by sequential 
projective measurements pT|. We call these "holonomic 
channels" . As we demonstrate below, the eigenvalues of 
the channel holonomy for smooth families of holonomic 
channels are directly related to the trace of the subspace 
holonomies, which reduces to Wilson loops in the sense 
of Ref. [1^ for closed paths. 



A. Unitary channels 

Here, we provide a detailed analysis of the holonomy 
for sequences of unitary channels, each of which acts on 
a Z?-dimensional state space. Unitary channels are in 
a sense the simplest type of channels, since their Kraus 
representations consist of a single unitary operator (K — 
1). Consequently, the Kraus freedom is encoded in a 
phase factor, and it follows that the holonomy in Eq. 
(|9]) is also a phase factor, which we denote 7ch- For a 
sequence £^^,£2 j ■ • • t^n '^^ unitary channels represented 
by the unitary operators C/i, C/2, ■ ■ • , C^jv, the holonomy 
7ch can be written 

7ch(fr,...,f]^) = <I>(Tr(f/^[/w))$(Tr([/;^[/Ar_i)) . . . 

x<i>(Tr(C/t;7i)), (38) 

where $(-z) = z/\z\ for any nonzero complex number 
z. Note that -fchi^i , ■ ■ ■ , ^n) gauge invariant in the 
sense that it is unchanged under the gauge transfor- 
mation Un e*""C/„ for arbitrary real numbers a„, 
n ^ 1,...,N. In particular, a gauge transformation 
Un — > Un such that 

HT^T{ul+,Un)) > (39) 



yields 



7ch(fi",...,f]^) = $(Tr(;7|[/jv)), (40) 



which is the channel holonomy in the parallel gauge (cf. 
Eq. (0) 

Consider now a smoothly parametrized family of uni- 
tary channels £f{p) — U{s)pW{s), s G [0,1], and con- 
sider the equation of motion iU{s) = H{s)U{s), H{s) 
being the "Hamiltonian" of the system {h — 1). This 
family defines the path C : [0,1] 3 s 1^ . To the first 
order in Ss we have TrlW {s + 6 s)U{s)] = D + idsTrH{s), 
and ^{Tt[W{s + Ss)U{s)]) = e-p*"'^''^(^' . In the 5s ^ 
limit, Eq. (|38p thus becomes 



(41) 



7ch(C) = $(Tr[C/t(G)C/(l)]) 



X exp 



D 



dsTrHis) 



which resembles the Aharonov-Anandan geometric phase 
Q with the "dynamical phase" —j^J^dsTTH{s) re- 
moved from the "total phase" arg$(Tr[f7'''(0){7(l)]). 
Note also that this dynamical phase coincides with the 
one in Eq. (14) of Ref. [23] for the maximally mixed 
internal state -^1 in an interferometer. 

The channel holonomy related to the family of unitary 
channels can alternatively be calculated using the paral- 
lel transport gauge, as developed in Sec. |TVl The parallel 
transport condition in Eq. (j25p reduces to the require- 

ment that Tr[J7 (s)J7(s)] should be a real number for all 
s. However, from the unitarity of U{s) it follows that 
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this number can only be purely imaginary. We thus find 
the parallel transport condition 



Tr[U {s)U{s)] = 0. 



(42) 



This condition is satisfied by the solution of the equa- 
tion of motion iU{s) = {H{s) - ^TT[H{s)]i}U{s). The 

unitaries U (s) and U (s) are related by the gauge trans- 
formation 



U{s) U{s) = U{s) exp 



DJo 



ds'TrH{s'). 



(43) 



The holonomy in the parallel transport gauge takes the 
form 



T,,h(C) = <i>(Tr[C/t(0)[/(l)]), 



(44) 



which can be seen by inserting Eq. into Eq. (|1T|) . 

Let us finally consider unitary channels for a single 
qubit (i.e., D = 2). We claim that one-qubit holonomies 
only can take the values ±1. To see this, we first note 
that U(2)=SU(2)xU(l). The U(l) parts of the sequence 
cannot contribute to the holonomy due to their cyclic 
appearance on the right-hand side of Eq. ((38)) . The claim 
then follows from $(Tr[SU(2)SU(2)]) = $(Tr[SU(2)]) = 
±1, since Tr[SU(2)] is a real number ^291]. 



B. Holonomic channels 



Consider a smoothly parametrized decomposition 
Hq — 0fcLi Tikis') of a D-dimensional Hilbert space Hq, 
for s' e [0, s]. Assume dim7ife(s') ~ Dk is constant for 
all k on the interval [0, s] and let Pk{s') be the projection 
operator onto 7ifc(s'). The quantities 



Tkis) = lim Pk{s)Pk{s- 

5s— i-0 



6s)...Pk{0) 



(45) 



can be expressed as [21| 

r,(s) = ^[Pe/o ^^(^')'^^'].,|af)(s))(af (0)1, (46) 

where P denotes path ordering, the anti-Hermitian ma- 
trix Ak{s') has elements [Ak{s')l,j = (df ^(s')|af ^(s')>, 

and {\a['^\s'))}f^i is an orthonormal basis of Tik{s'). A 
holonomic channel is defined as 



(47) 



The condition for trace preservation r[,(s)rfc(s) = i 
can be shown to be satisfied by using that Pe-^o ^'''^^ ^'^^ 
in Eq. I|46[) is a unitary matrix. 

Let us now examine the holonomy corresponding to a 
curve C : [0,1] 3 s £}}°^ of holonomic channels. The 
parallel transport condition in Eq. ([^S]) is satisfied if the 



matrix with elements Tr[rJ,(s)r;(s)] is Hermitian. By a 
direct use of Eq. ^ we find that Tr[f ^(s)r;(s)] = 
for all k,l and s. We can thus conclude that {rfe(s)}fc 
is parallel transported. Consequently, the holonomy is 
given by 

J7ch(C) = <f{To,i). (48) 
Here, Tg^i is a matrix with elements 
[To,i]ki = Tr[Ffc(0)rKl)] 

= c5.,^[Pe/o^^(^)'^1.,(af (0)|ap)(l)) 

= SkiTiiUgiCk)], (49) 

where U g{Ck) is the holonomy of the path Ck in 
the Grassmann manifold Q{D;Dk), i.e., the space of 
-Dfc-dimensional subspaces of an Z?-dimensional Hilbert 
space. Thus, the holonomy takes the form 

(7ch(C) = diag{<i>(Tr[J7g(Ci)]),$(Tr[J7<,(C2)]), 

...,<I>(Tr[i7g(C^)])}. (50) 

Note that if Cfc is a closed path, [To,i]fcfc equals the Wilson 

loopTrPe^'^.^'"'^^''^ 28]. 

A possible way to implement the holonomic channels 
is to approximate them by sequences of projective mea- 
surements. We discretize the interval [0, s] with step size 
Ss. We first perform the measurement {Pko{0)}ko fol- 
lowed by {Pki{Ss)}ki up to {Pfcjv(s)}fc„, where s = N6s. 
Discarding the outcomes of these measurements the re- 
sulting operation on the input density operator p reads 

P ^ PkAN6s)...PkMpPko{0)---PkANSs) 
= ^ Pfc (NSs) ...Pk {0)pPk (0) . . . Pfc {N6s) 

k 

+n{p) = M{p), (51) 

which is a channel. Here, 

np) = E PkANSs)...PkMp 

k() .. . .,kj\[ ^K, 

xPkM---PkANSs), (52) 

where /C is the set of values of fco,...,fcjv exclud- 
ing those where kg = ... = fcjv- Note that both 
Efc PkiNSs) . . . Pk{0)pPk{0) . . . PkiNSs) and 7^(p) corrc- 
spond to completely positive maps. In the Ss limit 
we have 



\im M{p)^£':,°Hp)+ lim n{p). 



(53) 



Now, A4 and are trace preserving from which it fol- 
lows that TZ must vanish in the Ss —* limit. Thus, we 
have 



limM{p)^£':,°Hp), 

(5s— *0 



(54) 



which concludes our demonstration that the holonomic 
channels can be approximated by sequences of projective 
measurements. 
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VI. PHYSICAL REALIZATION 

Given a physical system and a sequence of operations 
£i, . . . acting on this system, one might ask: what 
is the physical significance of the channel holonomy, and 
how should it be measured? Strictly speaking, these two 
questions have no well defined answer, given how we have 
constructed the channel holonomy. We have defined it as 
a change in Kraus representation, but the Kraus rep- 
resentation as such has no direct physical significance. 
Hence, there is no immediate way to attach measurable 
quantities to the channel holonomy. In other words, up 
to this point the channel holonomy has been a mathe- 
matical construction, rather than corresponding directly 
to a physical object or operation. However, we shall here 
obtain such a correspondence within the context of inter- 
ferometry. Another related question is, what is the object 
that the channel holonomy Uch represents? As Uch is a 
gauge-covariant unitary matrix, it seems intuitively rea- 
sonable that it should be a matrix representation of a 
gauge-invariant object. Our interferometric construction 
will provide precisely such an object. 



A. Yet another parallel transport 

Every channel £ on a Hilbert space TLq with dimension 
D can be obtained via a joint unitary evolution on the 
system and an ancillary system with Hilbert space Tia, 
according to 

E{p)^Txa[U{p(^\a){a\)U^], (55) 

where U is & unitary operator on Tiq ® Ha and \a) G Ha 
is normalized. If U is an arbitrary unitary operator on 
Ha, then both U and (1 (g) U)U give rise to the same 
channel £ [39| . If we regard the set of unitary operators 
U representing £ as a. fiber, the transformation Z^/ ^ (1® 
U)U can be seen as a gauge transformation resulting in 
a motion along the fiber. 

Consider now the sequence of channels £i, . . . ,£n, all 
with Kraus number K. We assume an ancillary space 
Ha of dimension K, and consider a sequence of unitary 
operators Ui, . . . ,Un on Hq (8) Ha, where Un represents 
£n via Eq. ((55|l . We regard the sequence of unitaries as 
parallel transported if 

TTq[Li„{i®\a){a\)L{l^^]>0. (56) 

Given an arbitrary sequence of unitaries Ui, . . . ,Un we 
can make it into a parallel transported sequence Un = 
{l(^Un)Un by choosing unitary operators Ui,...,Um such 
that 

Trq[Un{i®\a){a\)Ui^^] 

= U^Tvq[Un{i ® \a){a\)ui+^]Ul^^ > 0. (57) 



This requires that Trq[Z//„(i ® |a) (a|)W^^j^] is of rank K. 
It then follows that 

Un = Ui^{TTq[Uiii®\a)(a\)Ul])--- 

■■■<i>iTTq[UN-i{i®\a){a\)ul,]) (58) 

and Un = (1 <8) Un)Un- 

As may be seen, the above construction strongly re- 
sembles the channel holonomy in Sec. |TT1 as well as the 
Uhlmann holonomy in Sec. IIIII We shall see that this 
is not a mere coincidence, but that the above construc- 
tion enables us to obtain the channel holonomy within 
the context of interferometry. With this purpose in mind 
we review some of the concepts and tools that are useful 
for the analysis of quantum operations in single-particle 
interferometry. A more thorough account of these theo- 
retical tools can be found in Refs. [13, [3] • See also Refs. 
[31I m, [s^l for related material. 

B. Operations in interferometry 

Suppose that we have a single particle with some in- 
ternal degrees of freedom (e.g., spin or polarization) with 
Hilbert space Hq, and that the particle can propagate 
along two separated paths (e.g., the two paths in a Mach- 
Zehnder interferometer). These two paths correspond to 
two orthonormal vectors |0) and |1) spanning the "spa- 
tial" Hilbert space Hp = Sp{|0),|l)}. Thus, the total 
Hilbert space of the particle is Hp (S> Hq. 

Now, suppose that we have two operations Aq and Ai 
acting on the internal degrees of freedom of the particle. 
Let Aq operate on the particle if it passes path 0, and let 
Ai operate on the particle if it passes path I. The ques- 
tion is, what channels A acting on elements of C{Hp®Hq) 
would be compatible with the channels Aq and Ai act- 
ing in each path? (To be more precise, we require A to 
be a channel such that A(|0)(0| ® p) = |0)(0| ® Ao{p) 
and A(|l)(l| (g) p) = |1)(1| (g) Ai(p), for aU density opera- 
tors p on Hq.) The answer is as follows. Let {F„}„ and 
{W^m}m be linearly independent Kraus representations of 
the channels Aq and Ai , respectively. Then [l^ [ij] 

A(a) = |O)(O|0Ao((O|a|O)) + |l)(I|® Ai((l|a|l)) 
+ |0)(1|®^C„„„V;(0|<7|1)W^^ 

+ \l){0\®J2Cn,rnWm{l\a\0)V^ (59) 

nm 

where the matrix C satisfies CC^ < I and a S C{Hp (g 
Hq). We refer to A as a "gluing" of the two channels Aq 
and Ai , and to the matrix C as the "gluing matrix" . 

All gluings can be obtained using a shared ancillary 
system between the two paths Let V^^^ and V^^^ 

be unitary operators representing the channels Aq and 
Ai via Eq. (I55p . where we assume the same ancilla. On 
the combined system of the two paths, the system, and 
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the ancilla, we can construct the unitary operator 

C/tot = |0)(0|®V(°) + |l)(l|®V(i). (60) 

It turns out that 

A(a) - TTa{Utota®\a){a\Ul^) (61) 

is a gluing of the channels Aq and Ai, and moreover, that 
every gluing can be obtained by the appropriate choices 
ofV(o) and V(i) ^E^. Hence, although various choices 
of V^^-* and V*-^-* represent the same channels Aq and Ai, 
respectively, these choices may nevertheless result in dif- 
ferent gluings of Aq and Ai. 

C. Physical realization of the channel holonomy 

1. Interferometric parallel transport procedure 

We now consider the Mach-Zehnder setup. Let the 
particle start in path and internal state p (g) |a)(a|, af- 
ter which we apply a 50-50 beam splitter on the spatial 
degrees of freedom, acting as a Hadamard gate H on the 
spatial degrees of freedom regarded as a qubit. There- 
after, we apply JJtot in Eq- ((60)) on the total system fol- 
lowed by a variable unitary operator V^^^ on the ancillary 
Hilbert space in path and a variable unitary operator 
V^^"^ in path 1, i.e., the unitary operator 

F= |0)(0|®i®V^(°) + |l)(l|(»i®V^(^^ (62) 

Finally, we apply a second beam splitter and calculate 
the expectation value of the projector |0)(0| (g) 1 (g) 1, i.e., 
the probability p to find the particle in path 0, yielding 

p = Tr[(|0)(0|®i® i)H™tot 

xH(|0)(0|®/9® \a){a\)m^Ul^FW] 

= i + iRcTr{Wo)TrjV(°)(p'» |a>(a|)V(')V^'^^}- 

(63) 

If we regard the unitary operator V^^^ as fixed, then we 
find that the maximum probability p is obtained when 

1/(1) = y('')$(TrjV(°)(p® |a)(a|)V(i)^]). In the special 
case where p = jjl, i.e., the initial internal state is max- 
imally mixed, we thus find that 

yW = -i/(o)$(Tr,[v(°)(i (g |a)(a|)V(i)^]) (64) 

maximizes the detection probability p. 

Comparing with Eq. ([Ff|) we find that the parallel 
transport procedure can be implemented using this in- 
terferometric approach. We have a sequence of unitaries 
Hi,... Mn- We let V(°) = Ui and V^^) = U2, and choose 
the initial unitary operator V^") = Ui arbitrarily. Next, 
we let the input internal state to the interferometer be 
the maximally mixed state p = j^l and vary V^^^ until 



maximal detection probability is obtained. The maxi- 
mum is achieved when V^^'' = U2. In the next step we 
repeat the process but let V^"^ = U2 and V'^^ = U3 and 
we let F'-"-' = U2, i.e., the unitary operator we obtained 
in the first step. We then vary F^^) until we obtain max- 
imal detection probability. We repeat this procedure up 
to Un- 

One may note that the parallel transport procedure 
just described is purely operational, in the sense that it 
is achieved as a result of optimizing detection probabili- 
ties. This is analogous to the approach to the Uhlmann 
holonomy as developed in Ref. [33|. 

2. Realizing the channel holonomy Uch as a gluing between 
the end point channels 

After the A^th step of the parallel transport procedure, 
we can construct the channel holonomy as a measurable 
object in the interferometer. We let V'^-* = Un, V^^"-* = 
f/AT, V(i) =Z^i,andt/(i) =Ui. If we discard the ancillary 
system the resulting channel on C{T-lp ® TLq) is 

Afinai(a) = |0)(0| ®£:jv((0|a|0)) + \l) {l\ £i{{l\a\l)) 

+ mi\®Q{cj) + \l){Q\®g^a), (65) 

where 

g{a) = Tr,[(i ® C/w)ZiAr((0|t7|l) ® \a){a\)U\{l ® Ul)]. 

(66) 

We shall now see that the channel Afinai is directly related 
to the channel holonomy. 

Consider a unitary operator lAn representing via 
Eq. ([55)1 . Choose an arbitrary but fixed orthonormal 
basis {\ak)}k=i of for the evaluation of Tra. For such 
a choice, the operators E]^ = {ak\Un\a) form a linearly 
independent Kraus representation of the channel £„. We 
find 

g{a) = Y^{ak\ulUM\au')E^,{Q\a\l)El\ (67) 
fc.fc' 

One can furthermore show 

(afc|Tr,[Z^„(i® |a)(a|)Z^t+J| (68) 

(Note the reordering of k and k' between the left- and 
right-hand side.) As can be seen from the above equation, 
Ti q^n{^®\a){a\)ul^^i] is of rank K if and only if r„+i,„ 
is of rank K . Equation implies 

(afe|$(Tr,[Z^„(i ® \a){a\)Ui^^])\ak.) 

= [$(r„+i,„)]fc.fe, (69) 

which can be combined with Eq. (j58p to give 

{ak\UlUN\ak') = [<^{TN,N-i).-.HT2,i)]k'k 

= Zu'k- (70) 
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It follows that 
Afinai(a) = |O)(O|0£:jv((O|a|O)) + |l)(l|®5i((l|a|l)) 

K 

+ |0)(1|® Zu>^E^,{Q\a\l)El^ 

k,k' = l 
K 

+ |1)(0|® El,{l\a\Q)Efzl^,, (71) 

k,k' = l 

where the two Kraus representations and {E^^^ 

are free and independent of each other. We may there- 

TV 

fore choose a Kraus representation {£'; }; of the chan- 
nel £m, such that the overlap matrix with elements 

Tr{Ef Ej,) is positive definite, i.e., when the Kraus rep- 
resentation of £jv is parallel with the initial Kraus repre- 
sentation, which we know is possible if the overlap matrix 
[Ti^N]ik = Tt{EI''E^) is invertible. We thus substitute 
Ek' = Efe^f MTi,N)]kk' into Eq. (ED), which yields 

Afinal(a) - |0)(0|®£jv((0k|0)) + |l)(l|®5i((l|a|l)) 

K 

+ |0)(1|® [Uci.]k'kE^,{0\a\l)El^ 

k,k' = l 

+ |1)(0|® Y ^.Mlk|0)< [Ul^]kk'. 

k,k' = l 

(72) 

Thus, the channel holonomy C/ch is nothing but the glu- 
ing matrix describing the gluing of the two end point 
channels £i and £iq, with respect to the choice of paral- 
lel Kraus representations of these two channels. It is to 
be noted that a gauge transformation E], ^lEfVik 

iV N 

imply the same transformation Ej^ ~^ ^lEi Vik due 
to the assumption of parallelity between the end points. 
One may also note that the gauge covariance of C/ch, as 
described in Eq. is necessary for Afinai to be gauge 
invariant. Another way to put this is to say that U ch 
in some sense is the matrix representation of the gluing 
Afinai, and as a matrix representation of a gauge- invariant 
object we expect C/ch to be gauge covariant. We also 
point out that in the special case where the two end point 

channels coincide, i.e., Sn = £i, then Ef. = El. 

The above construction of the channel holonomy as a 
matrix representation with respect to two parallel Kraus 
representations is analogous to the construction of the 
open-path holonomy in Ref. (2l| . This is perhaps best 
seen in Eq. (12) in Ref. [2l[, where the open-path holon- 
omy C/g can be expressed asF = Y^kii'^ 9]ki\ak{l)) {ai{0)\. 
Here {|aA;(l))}fc is a frame which is parallel to the frame 
{|afc(0))}fc in a sense that resembles the parallelity of 
Kraus representations. (Compare, e.g., Eq. ([3]) of the 
present paper with Eq. (8) in Ref. |2ll|.) 

We finally note that the gluing matrix is in princi- 
ple possible to measure using interferometric setups . 



Hence, the channel holonomy resulting from the paral- 
lel transport procedure in Sec. IVI C II is a measurable 
object. 

VII. THE CASE OF SMOOTH 
PARAMETRIZATION 

In Sec. IIVI we considered the smooth version of the 
iterative parallel transport in Sec. HIl Here we perform 
the analogous transition for the ancillary construction in 
Sec. IVI Al finding conditions for parallel transport. 

Consider a smooth family of unitary operators U (s) 
acting on the combined system and ancilla Ti.qiS)'Ha, thus 
generating a smooth family of channels 

f,(p) = TTa[Uis)p® \a){a\U^{s)]. (73) 

When can we say that the family U (s) is parallel trans- 
ported? By recalling Eq. ([56]) it appears reasonable to 
require that 

Trg[Li{s){i(g)\a)(a\)U''{s + 5s)] > (74) 

to the first order in Ss in the limit of small Ss. We find 

TTg[Uis){i(S)\a){a\)U''{s + 6s)] = Q{s) + 6sTZ{s), (75) 

where 

Q{s) = Trg[U{s){i(g>\a){a\)U^{s)], (76) 

7^(.s)- TYg[U{s){i(g>\a){a\)u\s)]. (77) 

Note that Q{s) > 0. If we further assume that 

Q{s) > 0, (78) 

then we find that the parallel transport condition is sat- 
isfied when TZ{s) is Hermitian, i.e., 

TTq[U{s){i<»\a){a\)u\s)] 
^Tvg[U{s)il(g,\a){a\)U\s)]. (79) 

Suppose we have a smooth family of unitaries U (s) 
that is not parallel transported. The question is, under 
what conditions we can make it parallel transported by 
multiplying with a unitary on the ancillary space 

U{s) ^[i(g>U{s)]U{s), (80) 

i.e., a gauge transformation that leaves the path of chan- 
nels £s in Eq. ([73]) unchanged? Here we elucidate when 
it is possible to find a time-dependent Hamiltonian H(s) 
on Tia that generates U (s) via a Schrodinger equation. If 
we substitute Eq. ((80)l into Eq. (|79p we can rewrite the 
result as 

TZ{s) - TZHs) = Qis)U\s)U{s) + UHs)U{s)Q{s). (81) 
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If we now consider an anti-Hermitian operator A{s) gen- 
erating U{s) via U{s) = U{s)A{s), and substitute into 
Eq. dm) we find 7^ - 7^t = AQ + QA. Similarly as in 
Sec. IIVI we can use Theorem VII. 2. 3 in 24] to conclude 
that there exists a unique anti-Hermitian operator A{s) 
solving this equation. Hence, we can conclude that under 
the assumption Q(s) > 0, there exists a unitary family 
U{s) creating a parallel transported family U{s) via Eq. 

(EOl). 

It is to be noted that we have described a kind of 
two-step procedure. First, the system and ancilla evolve 
jointly according to U{s), which is not parallel trans- 
ported. In the second step, we modify the state by uni- 
tarily evolve the ancilla according to U{s). Hence, we 
have so far not obtained a joint Hamiltonian on the sys- 
tem and ancilla that generates the parallel transported 
family U{s). However, the latter can in principle be ob- 
tained since the Hamiltonian H{s) = ilAU generates U 
via the Schrodinger equation. Hence, the parallel trans- 
ported evolution can in principle be tailored through a 
combined evolution on the system and ancilla. Finally, 
we note that the parallel transport is not obtained oper- 
ationally in this smooth case, in the sense of the discrete 
case in Sec. IVI CI Whether such a "interferometric par- 
allel transport procedure" is possible in the smooth case 
we leave as an open question. 



VIII. THE 477 EXPERIMENTS 

Bernstein [34*1 and Aharonov and Susskind 'ss'l pointed 
out the possibility to observe the 47r spinor symmetry of 
a spin-i particle. The essence of their argument was 
that this symmetry may show up in the relative phase, 
say between two paths in an interferometer, one in which 
the spinor is rotated and one in which it is kept fixed. 
Subsequent neutron interferometer experiments [isl . [iGj 
were carried out to confirm this prediction. 

These experiments used unpolarized neutrons that 
were sent through a two-beam interferometer, in which 
one beam was exposed to a time-independent magnetic 
field B in the z-direction, i.e., B = Bi. In the weak-field 
limit, one can show that this corresponds to the family 
of unitary operations U{s) = e^^^"' , s e [0, tp], where 
ip (X B. A An periodic ip oscillation in the intensity in 
one of the output beams was observed by varying B. 

Let us consider the standard interpretation of these 
experiments. The spinor part of unpolarized neutrons 
is described by the density operator ^1, which can be 
decomposed into an equal-weight mixture of any pair 
of orthogonal pure spin-i states. Let the orthogonal 
vectors and represent such a choice of states. 

These vectors evolve into U{(p)\ip) and U {ip)\'ip^) , re- 
spectively, under the action of the magnetic field. It fol- 
lows that for (fi — n27r, we obtain \tp) (— l)"!";/;) and 
\ip'^) (— l)"!^-*-). Thus, irrespective how we choose I V') 
and IV'"'"), they have the desired 47r periodicity needed to 



explain the experiments. 

Here, we put forward another interpretation of these 
experiments, based on the channel holonomy. We demon- 
strate that the 47r periodic oscillations seen in the experi- 
ments can be interpreted as a An periodicity of the chan- 
nel holonomy in this case. We also show that the channel 
holonomy itself can be measured, by a slight modification 
of the setup in Refs. ,,15.. .16]. 

Let us first compute the holonomy of the path C asso- 
ciated with the family of unitary channels f " represented 
by the unitary operators U{s) = e~'^^'^' , s G [0,<p]. We 
obtain 

Tr[f7t(s)t/(s)]=-^TrM^0, (82) 

i.e., U{s) satisfies the parallel transport condition in Eq. 
dH]). Thus, we may use Eq. (gj) to deduce that 

7ch[q = $(Tr[[/(^)]), (83) 

since C/(0) = i. Explicitly, $(Tr[C/(^)]) = *(cosf), 
which is 47r periodic in if. 

Next, we analyze the experimental setup in Refs. 
[Tsl [T6j from the channel holonomy perspective. Let 
|0), |1) represent the two beam directions. These vectors 
constitute a basis for the spatial Hilbert space Tip. The 
experiment uses the standard interferometric sequence 
with an initial beam splitter, followed by an operation 
in the two paths, and finally a second beam splitter and 
a measurement. For the moment we focus only on the 
operation occurring in between the two beam splitters, 
but return below to the interferometer as a whole. Let 
cr be some arbitrary total state on Tip ® Hq . Apply to 
the beam the parallel transporting unitary family U{s), 
s e [0,(^], that represents the family of unitary channels 
along C. This results in the map 

A"(fT) = |0)(0|®£;((0|a|0)) + |l)(l|®£o"((lHl)) 
+ |0)(l|®7ch[C]C7(^)(0|a|l) 
+ |l)(0|®(l|a|0)C7^^)7ch[q*, (84) 

where we have used Eq. ([83]), and where U{ip) = 
^{TT[U^ip)])U{(p) is parallel to the initial unitary 
U{0) = i, i.e., Tr[U{ip)] > 0. We can identify this with 
Eq. ([7^ . Hence, A" is the gauge-invariant gluing be- 
tween the two unitary end point channels 6^—2 and 
Furthermore, the phase factor jch{C) is the corre- 
sponding gluing matrix with respect to the parallel Kraus 
representations C/(0) = 1 and U{(p) of the end point chan- 
nels. 

In the experiment in Refs. [l^, [l^ the input state to 
the interferometer was |0) (0|®^1 (unpolarized neutrons). 
If we assume that the beam splitters can be represented 
by Hadamard gates, the passage through the first beam 
splitter results in the state a = i(|0)-h|l))((0|-h(l|)«)ii. 
After the application of the channel A" on this state the 
particle passes through a second beam splitter, and the 
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probability to find the neutron in the beam, say, is 

Piv) = ^(l+ |Tr[C/(^)]|cos[arg7ch(C)]), (85) 

where we have used that Tr[C/(93)] = |Tr[J7((/3)] | . We 
see that pif) has period 47r in tp. Since the visibihty 
factor |Tr/7((/3)| = |cos^| is 2tt periodic in (p, the ob- 
served periodicity of the interference oscillations must 
originate solely from the An periodicity of the channel 
holonqm y 7n | i [C] . This concludes our demonstration that 
Refs. UMUM can be viewed as experimental realizations 
of a channel holonomy and its An periodicity. 

Finally, we wish to point out that the holonomy for 
this family of unitary channels could also be measured. 
It requires the following modification of the setup in Refs. 
[l^[la|: add a U(l) shift e'-^ to the 1 beam and maximize 
the output detection probability p by varying %. A direct 
calculation yields 

Pix) = ^(l + \^4UM\cos[x - arg7eh(C)]) , (86) 
which is maximal when x = arg7ch(C'). 

IX. CONCLUSIONS 

Investigations into quantum holonomy have yielded a 
unifying understanding of the geometry of some basic 
structures of quantum systems. These efforts have con- 
cerned not only the geometry of quantum states them- 
selves, but also how the twisting of subspaces, realized, 
e.g., as eigenspaces of some adiabatically varying Hamil- 
tonian, can be used to manipulate quantum states in a 
robust manner. 

In this paper, we have extended the notion of holonomy 
to sequences of quantum channels. The proposed quan- 
tity transforms as a holonomy under change of Kraus 
representations of the channels. We have shown that the 
channel holonomy is related to the Uhlmann holonomy 
for sequences of density operators constructed from 
the Jamiolkowski isomorphism 12]. We have delineated 
parallel transport and concomitant gauge potential for 
smooth families of channels. 

In addition to the relation to the Uhlmann holonomy, 
we have found some results that connect to other known 
quantum holonomies. For smooth sequences of unitary 
channels Jhe channel holonomy reduces to the phase shift 
in Rcf. [27] for unpolarized particles in an interferometer. 
Furthermore, we have analyzed a class of channels that 
have a direct relation to the subspace holonomies in Ref. 
PH . For smooth families of such "holonomic" channels, 
the channel holonomy is completely determined by the 
trace of the holonomies associated with paths in the space 
of subspaces (i.e., the Wilson loops if these paths are 
closed). 

We have demonstrated a physical realization of the 
channel holonomy in an interferometric setting, based on 



the idea of "gluings" of channels [l^, [13] • The realization 
consists of a gauge invariant object related to the chan- 
nel holonomy and a prescription for how this object can 
be used to extract the channel holonomy experimentally. 
Using this idea, we have been able to demonstrate that 
the neutron interferometer tests in Refs. [l^ [iB] of the 
An spinor symmetry can alternatively be interpreted in 
terms of the holonomy for unitary channels. Thus, one 
may be tempted to say that a particular realization of 
the channel holonomy already exists. 

To avoid technical complications, we have consistently 
made simplifying assumptions about the nature of the 
channels. We have restricted the analysis to sequences 
of channels with fixed number of linearly independent 
Kraus operators and we have assumed that all relevant 
matrices and operators have a well-defined inverse. If 
these assumptions are relaxed, one could consider ana- 
logues to the admissible sequences of density operator as 
considered by Uhlman n tsj] , or to the partial holonomies 
described in Refs. [2l|, [Sq . 

Let us end by a remark concerning the potential rel- 
evance of this work to holonomic quantum computation 
[lol . [ll[ . The key point with this type of quantum com- 
putation is that it is believed to be resilient to certain er- 
rors, such as those induced by open-system effects. Thus, 
in order to examine the resilience of holonomic quantum 
computation, it becomes important to have a useful no- 
tion of geometric phase or holonomy for open quantum 
systems. This has been addressed from different perspec- 
tives in several recent papers j59, 40, 41, 42, 4^^. 
Since open-system evolution may be described by quan- 
tum channels, it seems reasonable that the proposed 
channel holonomy, or some generalization of it allow- 
ing for variable Kraus number, might, in some way or 
another, play a role in the analysis of the resilience of 
holonomic quantum computation to open-system effects. 
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APPENDIX 

Consider channels £ and J- both with Kraus number 
K. Any choice of Kraus operators for this pair of chan- 
nels span /^-dimensional subspaces Cs and respec- 
tively, of the D'^-dimensional space of linear operators 
acting on the /^-dimensional state space. We show that 
a necessary and sufficient criterion for any T correspond- 
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ing to £ and T to have rank less than K is that Cs and 
Ljr are not fully overlapping. 

Assume Cg and are partially overlapping. Then, 
one can choose a Kraus representation {Fk}k for T so 
that there exists a -Ffe' in this set lying in the orthogonal 
complement to Cg. It follows that Tk'i = Tt:{fI,Ei) = 0, 
VZ, and thus the rank of T is less than K . 

Conversely, if the rank of T is less than then T 
has at least one singular value that is zero. Consider the 
singular value decomposition T = UDV {D diagonal 
and [/, V unitary). Assume Dk'W = and consider the 
transformation Fk Fk = X)/ ■ This results in 

T = DV, which implies fk'i = Ti{fIEi) = 0, VL 
Thus, Fk' lies in the orthogonal complement to Cg. 

As an illustration, consider the following representa- 
tions 



{Eo,Ei} 
{Fo,F,} 
(Go,Gi) 



(87) 



of the phase flip {£) , bit flip (jF) , and amplitude damping 
(Q) channels for a qubit {D = 2). Here, the ct's are 
the standard Pauli operators (with (7_|_ — cFx + io'y) and 
Pe.Pf.Pg e [0,1]. Clearly, K{£) = K{T) = K{g) = 
K ~ 2. By inspection, we see that Fi and Gi both lie 
in the orthogonal complement to Cs. Thus, both Cj^ 
and Cg overlap partially with Cs and the corresponding 
T matrices have rank less than K = 2. Explicitly, one 
obtains 



Tf,e 
Te,g 



2V(1-Pe)(l-P/) 





/^/T^(l + 



Clearly, the rank oiTp^E is 1 ii Pe,Pf ^ 1 and (Cs and 
Cyr orthogonal) if Pe = 1 or pf = 1. The rank of Tg,e 
is if Pe = 1 and pg = 0, and 1 otherwise. On the other 
hand, 



G,F 



yr^(i + yr^ 







VPgV^ ~ Pf 

which has rank K = 2 if pf ^ 1 and pg ^ 0. 



,(89) 
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